Non-Markovian effects on the non-locality of a qubit-oscillator system 
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Non-Markovian evolutions are responsible for a wide variety of physically interesting effects. 
Here, we study non-locality of the non-classical state of a system consisting of a qubit and an 
oscillator exposed to the effects of non-Markovian evolutions. We find that the different facets of 
non-Markovianity affect non-locality in different and non-obvious ways, ranging from pronounced 
insensitivity of the Bell function to quite a spectacular evidence of information kick-back. 
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In 1935, E. Schrodinger formulated a thought exper- 
iment addressing some paradoxical implications of the 
Copenhagen interpretation of quantum mechanics when 
pushed to the realm of everyday experience [l|. By de- 
scribing a situation where the degrees of freedom of a 
'large' object are correlated in a quantum mechanical way 
to a 'small' quantum system, the paradox by Schrodinger 
(commonly referred to as the 'cat paradox') embodies a 
genuine example of the possibility to enforce quantum 
features beyond the microscopic domain. Notwithstand- 
ing its almost octogenarian history, the cat paradox still 
defies a full understanding of its implications [2( . 

The steady-pace experimental progresses in quantum 
control achieved in the last twenty years have been able 
to produce instances very close to the original formu- 
lation by Schrodinger and are expected to help signifi- 
cantly in the grasping of fundamental concepts such as 
the quantum-to-classical transition, as well as the devel- 
opment of quantum technological applications @ . States 
having the form 



\^) = (n,D) s0 + \i,-D) s0 )/V2, 



(1) 



where {|t) s , \l) s } are the energy eigenstates of a spin-1/2 
particle (a qubit) and \±D) Q are opposite-phase coher- 
ent states of a harmonic oscillator Q, are faithful in- 
stances of the situations envisaged in Refs. [l| and have 
been demonstrated in trapped-ion settings Hl|. They 
are accessible (or close to be such) in other experimental 
contexts involving the effective interaction between spin- 
like systems and mechanical oscillators [oJ3] or the all- 
optical generation of micro-macro states [8[ . In the first 
instance, one would consider effective two-level systems 
(such as neutral or artificial atoms) embedded in cavi- 
ties endowed with movable light mirrors (embodying the 
continuous variable (CV) subsystem). In the second one, 
the spin and CV parts are provided by different degrees 
of freedom of two distinct photonic information carriers. 
Both settings are able to engineer states having the form 
of Eq. (TTJ) and both allow for the reconstruction of the 
Wigner function of the CV subsystem. As it will be seen 
in the next Section, such ability is crucial to the assess- 
ment of the Bell test at the core of our investigation. Re- 



markably, the multifaceted interests in studying quantum 
superposition states analogous to Eq. ([T]) extend up to 
the assessment of environment-induced dynamical effects 
and their implications for the settlement, manipulation 
and protection of general quantum correlations. This is 
even more relevant when non-trivial environmental in- 
fluences of a non-Markovian nature, such as those due 
to lack of divisibility dynamics and/or memory-keeping 
and feedback-inducing system-environment mechanisms 
are considered (9J. The working principles of such pro- 
cesses are still largely unexplored and are expected to be 
relevant in condensed matter set-ups involving artificial 
spins and mechanical modes. The experimental handi- 
ness of such states and the possibility to mimic the effects 
of non-trivial, memory-keeping environments in fully con- 
trollable linear-optics test-beds [l(| , make up for the pos- 
sibility to acquire knowledge on the true behavior of the 
quantum features of state (fTJ , when exposed to physical 
non-Markovian dynamics. 

Motivated by these arguments, in this paper we ad- 
dress the influences that non-Markovian dynamics giving 
rise to non-divisible maps have on the non-local nature 
of Eq. JT]) by studying two different configurations. First, 
we analyze the effects imparted by a spring-like coupling 
between the CV part of our state and an ensemble of 
quantum harmonic oscillators modelling quantum Brow- 
nian motion [llj . We then move to an effective post- 
Markovian dynamics of the spin-part only, as modelled 
by the master equation (ME) proposed in Re f. 113 and 
analyzed, for a single-qubit problem in Ref. 13]. Evi- 
dent signatures of non-Markovianity have been found in 
the trend of entanglement and discord [3] for an ini- 
tially quantum-correlated state of two harmonic oscilla- 
tors [15| . Here, not only we analyze a different figure of 
merit and form of quantum correlations, but we also ad- 
dress a radically different clas of states. We find that the 
behavior of the Bell function associated with a Brown- 
ian motion-affected superposition state of a qubit and an 
oscillator shows quite subtle features. First, Brownian 
motion affects the non-local nature of such state in quite 
a significant way when the cut-off frequency of the Brow- 
nian bath is much smaller than the natural oscillation 
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frequency of the CV subsystem, i.e. in the regime that 
would correspond to a strong non-Markovian limit: large- 
amplitude revival peaks are found, showing the kick-back 
mechanism that the memory-keeping environment can 
exert over the system. Second, the post-Markovian ME 
turns out to be unable to induce a non-monotonic decay 
of the Bell function. Yet, such dynamics is nondivisible, 
as it is straightforward to check, and as such it deviates 
from the prescriptions commonly accepted for Marko- 
vianity. Indeed, it is revealed as fully non-Markovian 
by the measure recently proposed by Rivas et al. [l6j . 
Remarkably, our study provides indirect evidence that 
the evolution of nonlocality in a quantum superposition 
state of a qubit and an oscillator is qualitatively similar to 
what would arise from the measure proposed by Breuer et 
al. [17[ , which is designed to point towards the back- flow 
of information from the environment to the system. 

The remainder of this manuscript is organized as fol- 
lows: In Sec. Uwe describe the formal tools that will be 
used in the core part of our analysis and discuss, very 
briefly, the case of Markovian evolutions. This will be 
used as a milestone for the comparisons with the explic- 
itly non-Markovian cases. In Sec. HH of the other hand, 
we present the key part of our study and address non- 
locality in a quantum superposition state of a qubit and 
an oscillator under non-Markovian dynamical conditions. 
Finally, Sec. IIIII is for our conclusions and outlook. We 
delegate the most technical parts of our work to two ap- 
pendices. 

I. TOOLS AND MARKOVIAN BENCHMARKS 

We start our study by describing the formal approach 
to non-locality that will be used throughout this work, 
which is similar to the one proposed in Ref. [l8| and used 
by Spagnolo et al. in Ref. [8]. In our thought experiment, 
the spin of the discrete- variable component of the system 
is probed along a direction n=(sin 9, 0, cos 9) of the Bloch 
sphere by the bi-dimensional operator 

<t(8) = sin 9 a x + cos 9 a z (2) 

with (T, (j—x,y,z) the j-Pauli operator. As discussed in 
Ref. [19( , non- locality of the state of a CV system can be 
tested in the phase-space by using the dichotomic parity 
operator 

oo 

n = (-l) ft = X)(l 2ri ) ( 2 "l - \ 2n + !> ( 2n + !|) (3) 

n=Q 

acted upon by the displacement T>(/3) = exp[/3a 1 ' — (3* a] 
(/J £ C) | so as to form the displaced parity opera- 
tor f[(/3) = 2?(/?)IIX' t (/3). Here, h = a^a is the bosonic 
operator number of the harmonic oscillator whose anni- 
hilation (creation) operators is a (aJ) and \n) is a Fock 
state with n excitations. The key point of such phase- 
space approach is that (FI(/3)) = (n/2)W((3), where 



W((3) is the Wigner function associated with the state 
over which the expectation value of fl(/3) is calculated. 
Therefore, we can easily construct the correlation func- 
tion C(/3 ) 0)=(i/j\a(9)^fl(/3) \if>) from which we get the 
Bell-Clauser-Horne-Shimony-Holt (CHSH) function 

B(fi, 9; 0', e')=C(p', e')+C{P', 9)+C(p, 6')-C{p, 9). (4) 

Local realistic theories impose the bound |£>|<2, which is 
violated by quantum mechanics, although not maximally, 
when state is used and a judicious choice of the set 
of parameters {/?, 9; (3', 9'} is made [18]. This form of 
the Bell-CHSH test has been shown to be effective in 
revealing the non-local nature of correlations in a state 
such as \ip) [l8| (see also Spagnolo et al. in Ref. 0) 
and we thus believe it is a very appropriate tool for the 
understanding of the role that the interaction with an 
environment has on the non-local properties of such state. 

First, we set a benchmark by briefly addressing 
the case of a Markovian dynamics as encompassed by 
general amplitude damping (AD) and phase damping 
(PD)processes, both for the spin and the CV subsys- 
tem [20] - Regardless of the dimensionality of the sys- 
tem at hand, both cases are most effectively tackled by 
means of the operator-sum representation of a quan- 
tum channel. We call {^(t)} the set of Kraus oper- 
ators specifying the non-unitary process p=AD,PD, so 
that the evolved state at time t takes the general form 
p(t)= Y, k ^fe(*) W -^(t). The explicit Kraus oper- 
ators for AD and PD and the form of the correlation 
functions corresponding to the case where the spin or 
CV components are influenced by the environment are 
provided in Appendix A. Here, it is enough to study 
the behavior of the numerically optimized Bell function 
max{g ,3.0/ ^/} \B\ for the superposition state, which is 
shown in Fig. [1] against the probability V v of occurrence 
of the channel (introduced in Appendix A). While, ex- 
pectedly, affecting the CV with an AD channel results in 
a quicker decay of the Bell function as compared to the 
spin-affected scenario, phase damping seems to be obliv- 
ious to the dimensionality of the subsystem being influ- 
enced by the environment. This is likely to be due to the 
fact that, as D gets sufficiently large, (— D\D) ~ and 
the CV state is effectively encoded in the fictitious quasi- 
qubit {|±Z?)} and the PD channel in the spin and CV- 
affected scenario give very similar outcomes. Needless to 
say, as the AD mechanism depends critically on the num- 
ber of excitations in the parties entering a given state, 
such 'homogenization' does not take place in Fig. [T] (a). 
In both cases, however, a monotonic decay of max \B\ is 
observed (for easiness of notation we omit the variables 
over which the optimization is performed). 

II. NON-MARKOVIAN EFFECTS ON 
NON-LOCALITY 

Our aim is now to study non-Markovian dynamics, 
looking for evidence of information kick-back over the 
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FIG. 1: We study max \B\ under the effects of amplitude- and 
phase-damping mechanisms. In (a) [(b)] max|£> is affected 
by AD [PD] and studied against the probability Vao ("Ppd) 
for the environmental action to occur. In both panels, the 
solid curve (dot-dashed curve) is for the spin-affected (CV- 
affected) case. The horizontal line shows the local realistic 
bound and we have taken D = 2. 



evolution of the hybrid system at hand. We start by an- 
alyzing the CV-affected case and consider the harmonic 
oscillator (with frequency uo) as coupled to an V-mode 
bosonic environment according to the interaction Hamil- 
tonian -ffcv=%X^=i c j9 ® Qj> where Cj is the coupling 
rate between the j th mode and the CV subsystem and 
g is a dimensionless coupling strength and q (Qj) is the 
position-like quadrature of the harmonic oscillator (the 



j mode of the environment). This model is known to 
give rise to quantum Brownian motion [TlT ] . The effective 
evolution of the CV subsystem is thus regulated by the 
time-local ME 



dtPo=~ {i/K)Uo [a f a, po] +Am (po ) , 
where po is the density matrix of subsystem O and 



£ bm (-)=-A(t)[q,[q,-}}+3(t)[q,[p,- 
-ij(t)[q,{p,-}} + - 



(5) 



(6) 



that accounts for diffusion [at rates A(t) and 3(i)], damp- 
ing [at rate j(t)] and the renormalization of the frequency 
of the CV subsystem. We have introduced the bosonic 
momentum-like quadrature p of the O subsystem. The 
derivation of Eq. ([5]) does not require the rotating-wave 
approximation nor it invokes the Born-Markov one. Var- 
ious dynamical phases have been identified for entan- 
glement and quantum discord of a two-mode Gaussian 
states under such dynamics [l^. These include effects 
of entanglement sudden-death and revival, depending on 
the spectral/memory properties of the environment itself. 

Here, we focus on the case of weak coupling, Ohmic en- 
vironmental spectral density with a cutoff U) c , high tem- 
peratures and short time-scale limit. Under these condi- 
tions, the contribution of r(t) to the solution of the ME is 
negligible, and the elements of the vector of coefficients 
v(t) = (A(t) 3(f) r y(t)) T entering Cb m are given by the 
expressions [151 ] 



vj(t)-- 



2 2 

9 uqx 
~2(l + x 2 ) 



fi(T){ aj 



-bj sin 



(7) 



with T=uj c t, x=uj c /ujo, T the environmental tempera- 
ture, f(T)—(k B T /ftu c )(l 1 haj c /k B T) T [k B is the Boltz- 
mann constant], a=(ill) T and b—(~lxx) T . The 
Brownian ME can be exactly solved using a phase space 
approach [2l| that allows for the determination of the 
Wigner function W(j3) = J-[xt(z)], where J 7 ^] indicates 
the complex Fourier transform, z=(qp) T is the vector of 
quadratures and Xt(z) is the Weyl function associated 
with the initial CV state [4] . Explicitly 



Xt(z) = e 



X0 (e-^ 2 R-\t)z) 



(8) 



with xo(^)=Tr[e 4 (' 5 «-^) / 9 (0)] the Weyl func- 
tion of the CV state at t=0, T(t)=2 J* j(s)ds, 
R(t)=cos(LUot)Tt+ism(u)ot)<Ty and 



W(t)-- 



-r(t) 



-R(t) 



r{s '>R T (s)M(s)R(s)di 



with M(s) 



2A(a) -3(a) 
-3(s) 



R T (t), 
(9) 

Due to the weak coupling 



and high temperature assumptions, g should take small 
values, while k B T '/Huj c cannot be too small. For short 
time scales, we can set e ±r< -^~l [la ]. 

With this at hand, we determine the Wigner functions 
(/?) of the CV components of the density matrix as- 
sociated to the spin part p\ J (i,j=^,l). Fig. [5] shows 
max \B\ as a function of t in two different dynamical 
regimes. For a large number of environmental mode- 
frequencies (i.e. for a large cutoff oj c ) and a relatively 
large coupling strength, the numerically optimized Bell 
function shows a monotonically decreasing behavior that 
makes us lose evidences of non-local character of the su- 
perposition state very soon in time. Although by low- 
ering g we observe some very partial revival of max|;B| 
due to a kick-back of coherence into the spin-CV sys- 
tem, this is not sufficient to give rise to violate the local- 
realistic bound again. The trend changes dramatically 
for uj c <C wo and j<C 1, which bring us to the phase of 
"non-Markovian revivals": max|B| becomes a periodic 
function of time and shows slowly-fading peaks at which 
the hybrid Bell-CHSH inequality studied here is quite 
largely violated. Such oscillations, which are typical of 
non-Markovian dynamics, are related to the appearance 
of temporal regions where A(t) achieves negative values 
and are connected to the memory of the environmental 
system, which keeps track of the system's state and feeds 
this information back to it. 

The connection between non-Markovian revivals of 
non-locality and the memory properties of the environ- 
mental system appears to be reinforced by the study of 
a simple model for system-environment interaction. We 
now consider the spin-part s of the superposition state 
as coupled to a star-like collection of N s non-interacting 
spin-1/2 particles via the energy-preserving longitudinal 

Ising model H sp i n =HAJ2k=i <*z®®z,k [&z,k is the 2;-Pauli 
operator of the fc th environmental spin and A is the cor- 
responding coupling strength]. The evolution induced 
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FIG. 2: Numerically optimized \B\ plotted against r for D = 2, ksT/hujc — 25. In panel (a) we have taken x = 10 and g — 0.3 
(solid line), 0.1 (dot-dashed line) and 0.05 (dashed one). Panel (b): max|B| for x = 0.2 and g — 0.05. The straight line shows 
the bound imposed by local realistic theories. 



by H S pi n over s is exactly solvable due to its excitation- 
preserving, non-interacting nature. We assume the spin 
star as prepared in the maximally mixed state 1/2 ", 
which leads to the reduced dynamics of spin s [22| 



S) = p1Ho) lt) s (tl 



[cos (2? 



P l sH0) \i) s d\ 
)] N '(ptH0)ft) a (l\ + h.c. 



(10) 



with t s = At, where p l J(0) (i, j=t, I) are elements of den- 
sity matrix p s (0). The calculation of the spin-oscillator 
correlation function proceeds now along the lines shown 
above and leads to the analytic form 



C(6, /^)=e^ 2|,3| '{sin6'cos(4D/3 l )[cos(2T s )] 



cos0sinh(4D/3 r )}, 



(11) 



where /3j=Im[/3], and /3 r =Re[/3]. The oscillatory func- 
tion of time appearing in C(6,/3) immediately reveals 
the non-monotonic behavior that the corresponding 
Bell-CHSH function will exhibit, thus signalling non- 
Mar kovianity. Clearly, at T s =m/2 (rGZ), the decoher- 
ence factor [cos^Ts)]^ 3 becomes ineffective, regardless 
of the number of environmental spins, and the correla- 
tion function achieves the value corresponding to a pure 
qubit-oscillator superposition state. Therefore, the non- 
Markovian non-locality revivals are full and the only ef- 
fect of a growing size of the star is the narrowing of the 
revival peaks [23j . This results in shorter time windows 
where the Bell-CHSH inequality is violated. This analy- 
sis is displayed in the upper part of Fig. [3] (a) , where we 
assess the case of three different values of N s . 

An interesting remark is due, now, in relation to the 
aims of our study. This spin model for decoherence is 
known to provide a divergent value of the non-Markovian 
measure proposed by Breuer et al. in Ref. (l7j [which is 
unbounded]. This means that the dynamics experienced 
by s can never be described by a Markovian model. The 
model, in fact, is such that the trace distance of equa- 
torially antipodal states of s (upon which the measure 



in [17[ is built) is given exactly by | cos(2r s )| JVs . This is 
clearly shown in the lower part of Fig. [3] (a). Therefore, 
not only the non-Markovian revivals of non-locality per- 
sist in time due to the infinitely non-Markovian nature 
of the evolution, but their occurrence is clearly related to 
the changes in the trace distance, thus providing a clear 
connection between the kick-back of information from the 
star system to s and the revived non-local features of \ip). 

However, one should be careful in dealing with the re- 
lationship between non-Markovianity and non-locality, 
due to the multifaccted nature of the former statisti- 
cal phenomenon: it would be a mistake to identify an 
in-principle non-Markovian dynamics with the occur- 
rence of non-monotonic trends in max|S|. We illus- 
trate such point considering an important instance of 
post-Markovian dynamics proposed by Shabani and Li- 
dar in [12]. The model includes explicitly a memory ker- 
nel k(t) that renders the ME time non-local as 



d t p s (t) = L 



k{t')e Lt ' p s (t 



t')dt', 



(12) 



where L[-\ is the standard Markovian Liouvillian de- 
scribing dissipation at rate 70 induced by a thermal 
bath with mean occupation number n 0, H^J. We 
focus on the widely used form of the memory kernel 
k(t)=je ~ 7 *, a choice that guarantees complete positivity 
of the map for any value of (7, 70, n) and at any instant of 
time [H, [26| . The dynamics in Eq. (|12[) can be straight- 
forwardly solved using a rather technical approach dis- 
cussed in Refs. fUJI!?]]. For the sake of providing a self- 
contined presentation, we summarize the steps needed 
for such a derivation in Appendix B. We eventually get 
an analytic form of the density matrix of the evolved 
superposition state and, thus, the correlation function, 
whose expression is however too cumbersome to be re- 
ported here. Fig. [3] (b) shows the results of our quanti- 
tative study: we display the optimized Bell-CHSH func- 
tion against the dimensionless interaction time T s i—^^t 
and for a long/short environmental memory time as com- 
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FIG. 3: (a) Upper part: Time dependence of the numerically optimized Bell function for a quantum superposition state between 
a qubit and an oscillator with D—2 and whose s subsystem is affected by iV s =2,5, 100 non-interacting spins according to the 
model H spin- Here, r s =At is a dimensionless interaction time. Lower part: Behavior of the trace distance 8(T)=Tr\p s ,+ —p s ,- |/2 
(with p s ,±= |±} s (±| and a x |±) s = ± |±) s ) upon which the measure of non-Markovianity proposed in [17| is built. From top 
to bottom curve, we have N s = 2, 5, 100. (b) We show max \B\ against r s ; for 70/7=0. 05, 1, 10 (solid, dashed and dot-dashed 
line, respectively) and n = 0. The inset shows max|£>| against n for r s i — 1.6 and 70/7 = 10, 14.3 (solid and dashed line, 
respectively). The horizontal line shows the local realistic bound. We have D=2 in all plots. 



pared to the dissipation time Jq 1 - When 7^70, the 
environment has a long-time memory, thus pushing the 
dynamics away from Markovianity. Interestingly, this re- 
sults in larger temporal windows where the violation of 
the Bell inequality is observed. Viceversa, for a larger 
7, such window becomes shorter and the dynamics gets 
closer to a Markovian one. The resilience induced by 
a memory-keeping evolution turns out to be quite spec- 
tacular when studied against the thermal nature of the 
environment. As shown in the inset of Fig. [3] (b), while 
for 70/7=10 it is enough to have n=1.6341 to stop violat- 
ing the Bell-CHSH inequality, a small raise to 7o/7~14.3 
is enough to push the value of n at which max |£>|<2 to 
n~200. 

It is remarkable, though, that besides such induced 
robustness, no otherwise evident signature of non- 
Markovianity (such as the ripples highlighted previously 
and typically exhibited by other indicators of quantum- 
ness [151 ]') can be deduced from the analysis of the cor- 
relation function C(6,{3) associated with this case and 
of Fig. [3] (b). We believe that the reason for such a 
behavior, though, should be researched in the multi- 
faceted nature of non-Markovian dynamics. In fact, the 
Shabani-Lidar ME can be recast into the form of a time- 
dependent, time-local ME whose coefficients are not all 
positive, therefore signaling the break-down of divisibil- 
ity 28j and thus the impossibility to describe the evolu- 
tion in Markovian terms, according to the criterion put 
forward by Rivas et al. 16], although no kick-back of in- 
formation is possible (as witnessed by the fact that the 
measure proposed in |17| is strictly null, for such a dy- 
namical map). It is just intriguing that the Bell-CHSH 
function is able to experience the subtleties of such dif- 
ferences in such a striking way. 



III. CONCLUSIONS 

We have studied the behavior of the Bell-CHSH func- 
tion for a spin-oscillator superposition state under the 
influences of two models for non-Markovian dynamics, 
demonstrating the emergence of interesting features re- 
lated to the various facets with which non-Markovianity 
manifests itself. Care should be used in assessing them by 
means of diverse non-classicality indicators, as they may 
turn up to give mutually inconsistent evidences. We have 
shown that the Bell-CHSH function is particularly sensi- 
tive to the subtleties of a dynamics that, although does 
not appear to be characterized by a kick-back of informa- 
tion from the environment, is nevertheless not divisible 
and as such far from being Markovian. It will be inter- 
esting to address similar questions in models for open- 
system dynamics specific of condensed matter set-ups, 
such as the 1// noise that is expected to provide strong 
non-Markovian dynamical features. This will help us as- 
certain if the range of non-Markovian manifestations in 
non-locality tests is even richer than the predictions com- 
ing from our work. 
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APPENDIX A 



in this case reads 



Here we provide the explicit form of Kraus opera- 
tors for Markovian amplitude damping (AD) and phase 
damping (PD), and also the specific correlation functions 
corresponding to the case where the spin or CV compo- 
nents of our superposition state are affected by the envi- 
ronment. 



A. Amplitude damping 

The general Kraus operators for a d-level system ex- 
posed to an amplitude damping (AD) channel character- 
ized by the damping rate 7 are [291 ] 

■W*) = E V {n - k)/2 (l-v) k/2 \n-k)(n\ (A-l) 



with k=0, ..,d—l, Tj 



and [1 — rj\ k l 2 the probability 



that a Fock state \n) loses k excitations within time t. 



1. Spin damping 

For a spin- 1/2 particle we have d = 2 and the operator- 
sum decomposition consists of the following two elements 



M D {t) = 



v^ii)<ii, A? D (t) = yn^io><i|. 

(A-2) 



We call Pad^a/I ~~ V f ne probability that the system 
loses one particle up to time t and identify |0) (|1)) with 
state \l) (|t))- In this case, the correlation function for 
the Bell-CHSH non-locality test performed in the body of 
the manuscript takes the form (for simplicity, we assume 



2 277-1-e 8 ^ 



which is such that |C~f n (/3, 6, r?)|<l. 



cos 9+^/rje p s'md. 

(A-3) 



2. CV damping 

When considering a CV system, the summation in 
Eq. (|A-1I) extends to d — > 00. Assume that the system is 
prepared in a coherent state Such state is changed 
by the action of the fc th Kraus operator as 



Ai D (t)\0 = e 



\VvO, (A-4) 



showing that a coherent state remains such under the ac- 
tion of an amplitude damping channel, although its am- 
plitude is reduced. The Bell-CHSH correlation function 



(A-5) 



sin f 



B. Phase damping 

The general Kraus operators for a d-level system ex- 
posed to a phase damping (PD) channel characterized by 
the rate \i are [3C| 



n=0 



where r pd — fit is the rescaled interaction time. 

1. Spin damping 

For a spin-1/2 particle, the Kraus operators are 
A PD (^) = |0)(0|+e-^|l)(l|, 



(A-6) 



(A-7) 



^( Tpd)= Vl- e -^|l)(l|, 



where Vpr>= T f d is the probability that one exci- 

tation from the system is scattered by the environment. 
The Bell-CHSH correlation function in this case reads 



C S pm(/3j $7 Tpd)- 



-2(P+D) 2 



-(l-e 8 ^)cosfl+e- 2 ^-^sin< 

(A-8) 



2. CV damping 

Again, let us assume that the CV system is prepared in 
a coherent state Such state is changed by the action 
of the fc th Kraus operator as 

T k 00 ken 

A PD (r pd )\0 = ^£ e -i(«%W)ILf | n) . (A-9) 

v n— v 

In the spin-1/2 basis {|t)s, \i)s}, the Schrodigner cat 
state that has been affected only by a CV PD channel 
has the representation 



rev 



1 



N ++ ( Tpd ) N + _(r pd ) 
N-+(T pd ) N—( Tpd ) 



(A-10) 



with 



k—0 n,m— 

(nm) k (iD) n (jD*) 



\ri){m\, (i,j = ±). 

(A-ll) 
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Let us first calculate the expectation value of the projec- 
tion operator a (9) defined in the body of the manuscript 
over the state p^ v D . We have 



-, oo 2k oo 



with 



2 ^ k\ 

k—0 n,m— 



|2 } (nm) k D n {D*y 



(A-12) 



\Jn\m\ 

x {sin(9[(-l) m + (-1)"] + cos 6» [(-1) - 1]}. 

(A-13) 

The Bcll-CHSH correlation function is then found by cal- 
culating the expectation value of the displaced parity op- 
erator n(/3) as 



C™{p,9,T pd ) = Tr[n(/3)£(0,r pd )]. 



(A-14) 



During this calculation, we need to evaluate 

the expectation value of the parity operator 

(-l)"=^ =0 (|2n)(2n|-|2n+l)(2n+l|) over displaced 
Fock states as 



( m \D(f3)(-lfD\f3)\n), 
which can be readily evaluated using [3l| 



(A-15) 



(sMP)\r) = ^(Py- r e-^Ll s - r \\^), (s > r) 

(A-16) 

where L p l \x) is an associated Laguerre polynomial. We 
now take four different cases 

1. For m > 2n +1 with 2n > n. In this case we 
introduce the function 



Si 



-{li) m ' n e~^\-l) 2n '~ n [L ( Zr 2n '\\P\ 2 ) 



xi( 2 »'-«)(|^| 2 )+4--f- 1 )(|/3| 2 )4 2 »'+ 1 -™)(|^| 2 )]. 

(A-17) 



2. For 2n > (m,n), we define 



S 2 



V m\n\ 



2n'\ 



\2n —n. 



-\0\ 2 ]L (2n'- 
L m 



4. For n > 2n +1 and 2n > m, we introduce 

Si = y^(/3*) n - m e- | ' 9|2 (-l) 2n '- m [L I (2ll '- m) (|^| 2 ) 

x4^ 2 "' ) (i/?i 2 )+4r ,+i - m) (i/?i 2 )4r + 2 r" 1) (i/?i 2 )]- 

(A-20) 

With such definitions, the Bcll-CHSH correlation func- 
tion is finally obtained as 

Cr v D (l3 : 8,T pd ) = 

1 oo 2k oo I oo 2n' 2n' 2n / 

*E^f£ E E^ + EE^ 2 

oo oo oo 2n 

E E ns »- E E^ 

m— 2rc/ + l n— 2n' + l n— 2n' J r l m—0 

(A-21) 



APPENDIX B 

For the sake of providing a self-contained presentation, 
here we sketch a strategy to gather an analytic solution 
of the post-markovian master equation (ME) assessed in 
the body of the manuscript 12j. We follow the lines 
sketched in Refs. [13|,|27| and first introduce the operator 
eigenbasis {Qk} of the superoperator 



L{p s ) = 7o(« + l)(a-p s a + - -{a+a^,p s }) 
+ 7 n(o- + p s <5-_ - i{o-_<T + ,p s }) 



(B-l) 



with a± = (a x ± ia y )/2 and p s the general density ma- 
trix of the spin system. The elements of such eigenbasis, 
which is commonly referred to as the damping basis, are 
such that 



(B-2) 



with {Afc} the set of corresponding eigenvalues. It is 
straightforward to find that [27[ 



1 - 



(271+1) 



, Q 2 = <J Z , Q3 = ct+, Q4 = cr_ 



x4 2 "'-")(|/3| 2 )--^4 2 /+ 1 - m )(|/3| 2 )4 2 "'+ 1 -")(|/3| 2 )] 
In + 1 

(A-18) 



(B-3) 

with {Ai = 0, A 2 = 2A 3 ,4 = -27o(n+l/2)}. In the damp- 
ing basis, the density matrix of the system can be written 

as 



3. For (m, n) > 2n +1 we call 

* 3 = 4S=(/3) , ^ 2n '(/3*) n " 2n 'e- |/3|2 [4r 2 "' ) (i/3i 2 ) 



(B-4) 



fe=i 



yrrdn). 

(A-19) 



Substituting Eqs. (fB^2|) and (jB^4|) into the ME (10) in 
the main manuscript, the equation changes to 



£>i(t)-A(t)]&=o. 



(B-5) 
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where /3 k (t) = /„ k(t')X k e Xkt a k (t - t')dt' . From this, we 
easily extract the set of integro-differential equations 

a l (t)=/3 l (t) (1=1, ..,4), (B-6) 

which can be solved by resorting to Laplace transforms 
to get 



where X (s):=Lap[X (t)] is the Laplace transform of X(t) 
and Lap -1 stands for the inverse transform. The figures 
of merit studied in the main paper are then easily deter- 
mined from such solutions. 



ct p (t) = Lap 



1 



X p k(s — A p ) 



ap(0), 



(B-7) 
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